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We present a generic mechanism by which reproducing microorganisms, with a difTusivity that 
depends on the local population density, can form stable patterns. It is known that a decrease of 
swimming speed with density can promote separation into bulk phases of two coexisting densities; 
this is opposed by the logistic law for birth and death which allows only a single uniform density to 
be stable. The result of this contest is an arrested nonequilibrium phase separation in which dense 
droplets or rings become separated by less dense regions, with a characteristic steady-state length 
scale. Cell division mainly occurs in the dilute regions and cell death in the dense ones, with a 
continuous flux between these sustained by the diffusivity gradient. We formulate a mathematical 
model of this in a case involving run-and-tumble bacteria, and make connections with a wider class of 
mechanisms for density-dependent motility. No chemotaxis is assumed in the model, yet it predicts 
the formation of patterns strikingly similar to those believed to result from chemotactic behavior. 



Microbial and cellular colonies are among the simplest 
examples of self-assembly in living organisms. In na- 
ture, bacteria are often found in concentrated biofilms, 
mats or other colony types, which can grow into spectac- 
ular patterns visible under the microscope [TJ [5] . Also 
in the laboratory, bacteria such as E. coli and S. ty- 
phimirium form regular geometric patterns when they 
reproduce and grow on a Petri dish containing a gel such 
as agar. These patterns range from simple concentric 
rings to elaborate ordered or amorphous arrangements 
of dots [3H6]. Their formation results from collective be- 
haviour driven by interactions between the bacteria, such 
as chemotactic aggregation [6!, competition for food [S] 
or changes in phenotypes according to density |llj . The 
question as to whether general mechanisms lie behind 
this diversity of microscopic pathways to patterning re- 
mains open. 

Unlike the self-assembly of colloidal particles, pattern 
formation in motile microorganisms and other living mat- 
ter is typically driven by non-equilibrium rather than 
thermodynamic forces. Indeed, the dynamics of both di- 
lute and concentrated bacterial fluids is already known 
to be vastly different from that of a suspensions of Brow- 
nian particles. For instance, suspensions of active, self- 
propelled particles, have been predicted to exhibit gi- 
ant density fluctuations [TH [T3], which have been ob- 
served experimentally jl4| . Similarly, an initially uni- 
form suspension of self-propelled particles performing a 
"run-and-tumble" motion like E. coli has recently been 
shown theoretically to separate into a bacteria-rich and a 
bacteria-poor phase, provided that the swimming speed 
decreases sufficiently rapidly with density [15j . This is 
akin to what happens in the spinodal decomposition of 
binary immiscible fluids, but has no counterpart in a sys- 
tem of Brownian particles interacting solely by density- 
dependent diffusivity. (The latter obey the fluctuation- 
dissipation theorem, ensuring that the equilibrium state 



is diffusivity-independent.) Other non-equilibrium ef- 
fects, such as ratchet physics, have also been observed 
and used either to rectify the density of bacteria [l6l - fT8] 
or to extract work from bacterial assemblies [T^ . 

Some aspects of bacterial patterning show features 
common to other nonequilibrium systems, and a cru- 
cial task is to identify the key ingredients that control 
their development. In many equilibrium and nonequi- 
librium phase transitions an initial instability creates 
density inhomogcneitics; these coarsen, leading eventu- 
ally to macroscopic phase separation [5D]. The situa- 
tion observed in bacterial assemblies often differs from 
this; long-lived patterns emerge with fixed characteristic 
length scales, suggesting that any underlying phase sep- 
aration is somehow arrested. The strong diversity of bio- 
logical functions met in experiments has led to an equally 
diverse range of proposed phenomenological models [5]- 
7, 9 llj to account for such effects. Most of them rely 
on the coupling of bacteria with external fields (food, 
chemoattractant, stimulant, etc.), and many involve a 
large number of parameters due to the complexity of the 
specific situation of interest. The most common mech- 
anism used to explain the bacterial patterns is chemo- 
taxis [S]: the propensity of bacteria to swim up/down 
gradients of chemoattractants/repellants. This explana- 
tion is so well established in the literature for at least two 
organisms {E. coli and 5". typhimuridium [6 ) that obser- 
vation of similar patterns in other species might defensi- 
bly be taken as evidence for a chemotactic phenotype. 

Here we identify a very general mechanism that can 
lead to pattern formation in bacterial colonies and which 
may encompass a large class of experimental situations. 
This mechanism involves a density-dependent motility, 
giving rise to a phase separation which is then arrested, 
on a well defined characteristic length scale, by the birth 
and death dynamics of bacteria. For definiteness we 
will work this through for a particular model of bacte- 



2 



rial run-and-tumble motion, involving a swim speed that 
depends (via unspecified interactions) on local bacterial 
density. This gives pattern similar to those observed 
in experiments [51 However, the basic mechanism 
— density-dependent motility coupled to logistic popu- 
lation growth — is not limited to this example. Our work 
demonstrates that chemotaxis per se is not a prerequisite 
for observing what are sometimes colloquially referred to 
as 'chemotactic patterns'. 

It is indeed remarkable that density-dependent swim 
speed and logistic growth alone are sufficient to create 
some of the specific pattern types previously identified 
with specific chemotaxis mechanisms. In mechanistic 
terms, we find that the logistic growth dynamics effec- 
tively arrests a spinodal phase separation that is known 
to follow from a density-dependent swim speed [15] . Put 
differently, an initially uniform bacterial population with 
small fluctuations will aggregate into droplets, but these 
will not coarsen further once a characteristic length scale 
is achieved, at which aggregation and birth/death effects 
come into balance. Starting instead from a small inocu- 
lum, we predict formation of concentric rings which, un- 
der some conditions, at least partially break up into spots 
at late times 

To exemplify our generic mechanism we will start from 
a minimal model of run-and-tumble bacteria, that can 
run in straight lines with a swim speed v and randomly 
change direction at a constant tumbling rate fITl\'2'2\ . 
To this we add our two key ingredients: a local density- 
dependent motility, and the birth/death of bacteria, the 
latter accounted for through a logistic growth model. 
Of course, bacteria can interact locally in various ways, 
ranging from steric collisions |15| to chemical quorum- 
sensing [5] . (Indeed a nonspecific dependence of motility 
on bacterial density was previously argued to be central 
to bacterial patterning by Kawasaki [5].) Here we fo- 
cus on the net effect of all such interactions on the swim 
speed v{p), which we assume to decrease with density 
p. This dependence might include the local effect of a 
secreted chemoattractant (such as aspartate [3H5] which 
causes aggregation, effectively decreasing v) but does not 
assume one. 

In addition to their run-and-tumble motion, real bac- 
teria continuously reproduce, at a medium-dependent 
growth rate which ranges from about one reciprocal hour 
in favourable environments such as Luria broth to sev- 
eral orders of magnitude lower for 'minimal' media such 
as M9. In bacterial colonies patterns may evolve on 
timescales of days [6] , over which such population growth 
dynamics can be important. 

We now derive continuum equations for the local den- 
sity p{r, t) in a population of run-and-tumble bacteria, 
with swim speed v{p), growing at a rate a(l — p/po). 
The latter represents a sum of birth and death terms, in 
balance only at p = po. At large scales in a uniform sys- 
tem, the motion of individual bacteria is characterized by 
a diffusivity D{p) = v(p)'^T /d, where is the tumbling 
rate and d the dimensionality [211 [Hj. Crucially how- 



ever, a non-uniform swimming speed v{v) also results in 
a mean drift velocity V = —vtS/v |21j which here gives 
V = —D'{p)Vp/2 [IS]. This contribution is crucial to 
phase separation 15] and will again play a major role 
here. However this term has no counterpart in ordinary 
Brownian motion (even if particles have variable diffusiv- 
ity) and was accordingly overlooked in previous studies 
which relied on phenomenological equations involving a 
density-dependent diffusivity and no drift ^8i. 

Coupling the diffusion-drift equation for run-and- 
tumble bacteria, as derived in [151, with the logistic 
growth term, the full dynamics is then given by: 

= V-[I?e(p)Vp(r,i)] (1) 

+ap{r,t) (^1 - - «VV(r,0 

where the 'effective diffusivity' is 

V,{p)^D{p) + ^pD'{p) (2) 

This results from the summed effects of the true diffu- 
sive flux —D{p)\Ip and the non-linear drift flux pV . In 
Eq.([T]) we have also introduced a phenomenological sur- 
face tension k > 0, which controls gradients in the bacte- 
rial density. Such a contribution has been shown to arise 
when the speed of a bacterium depends on the average 
density in a small local region around it, rather than a 
strictly infinitesimal one [15 . Eq. ([T]) neglects noise, both 
in the run-and-tumble dynamics and in the birth/death 
process. The former noise source conserves density and 
should become irrelevant at the experimental time scale 
of days. On the other hand, the non-conservative noise in 
the birth and death dynamics may be more important, 
and we have verified that our results are robust to its 
introduction at small to moderate levels. Numerical sim- 
ulations of Eq. ([T]) have been performed with standard 
finite difference methods (although noise does require 
careful treatment, as in [23 HI]), with periodic bound- 
ary conditions used throughout. For definiteness, all our 
simulations have been carried out with v{p) — iigc"^''^^, 
where > is the swim speed of an isolated bacterium 
and A > controls the decay of velocity with density. 
The precise form of v{p) is however not crucial for the 
phenomenology presented here, and the instability anal- 
ysis offered below does not assume it. 

The logistic population dynamics alone would cause 
the bacterial density to evolve toward a uniform den- 
sity, p(r) = which constitutes a fixed point for the 
proposed model. Although this homogeneous configu- 
ration is stable in the absence of bacterial interactions, 
it has been shown [T5| that, without logistic growth, a 
density-dependent swim speed v{p) leads to phase sepa- 
ration via a spinodal instability whenever dv /dp < —v/p. 
By Eq. ^ this equates to the condition < 0, and it 
is indeed obvious that the diffusive part of Eq. ([ij is 
unstable for negative I?o- It is important, clearly, that 
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Dc can be negative although D is not. This holds for a 
much wider class of nonequilibrium models than the one 
studied here; we return to this point at the end of the 
paper. 

For the choice of v{p) made in our simulations, we 
have 2?c = D{p)[l — pX/2] and the flat profile will thus 
become unstable for po above 2/ A. We have confirmed 
this numerically, and find that upon increasing po, the 
uniform state becomes (linearly) unstable, evolving in a 
ID geometry into a series of "bands" of high bacterial 
density separated by low density regions. Depending on 
the parameters, this transition can be continuous (super- 
critical), with the onset of a harmonic profile whose am- 
plitude grows smoothly with po, or discontinuous (sub- 
critical) with strongly anharmonic profiles (see Fig. f ). 




Figure 1: Growth of the instability in the supercritical (left) 
and subcritical cases (right). The three lines correspond to 
three successive times. A small perturbation around po (red 
line) growth toward harmonic or anharmonic patterns in the 
supercritical or subcritical case, respectively. Left: Super- 
critical case {a — K, — 0.01, A — 0.02, po = 15, -Do ~ 1, times: 
10^, 10^, lO"*). Right: Subcritical case (« = «: = 0.005, 
A = 0.02, po = 11, Do = 1, times: 3. 10^ 3. 10^ 10^) 



The transition to pattern formation arising from 
Eq. is a fully nonequilibrium one: it is not possible to 
write down an effective thermodynamic free energy which 
would lead to this equation of motion. Nonetheless, it is 
possible to understand why the birth/death process effec- 
tively arrests the spinodal decomposition induced by the 
density-dependent swim speed. The latter tends to sep- 
arate the system into high and low density domains with 
densities either side of po- (Without the logistic term, 
these would coarsen with time.) Bacteria thus tend to 
be born in the low density regions and to die in the high 
density regions. To maintain a steady state, they have to 
travel from one to the other: balancing the birth/death 
terms by the diffusion-drift transport flux between the 
domains then sets a typical scale beyond which domain 
coarsening can no longer progress. Were any domain to 
become much larger, the density at its centre would soon 
regress towards po, re-triggering the spinodal instability 
locally. (This is closely reminiscent of what happens in a 
thermodynamic phase separation when the supersatura- 
tion is continuously ramped [25 .) 

To better understand the onset of the instability, let 
us linearize Eq. (|T]) around p(r) = po and work in Fourier 
space. Defining p(r) = po + J2q ^Pq exp(iq • r) yields: 




Figure 2: Three plots oi ^q{q) for \V^{po)\/ s/an = 1;2;3. At 
the transition, only one critical mode q = Qc is unstable. 



The flat profile p = po is thus stable if < for all 
q and is unstable otherwise. From the expression for 
Vcipo), Eq. ([2]), one sees that instability occurs if 



poi^'(po) 
- 2Dipo) - 



and -^>2 
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At the onset of the instability only one mode is unsta- 
ble, with wavevector Qc — ^2q;/|I?o(po)|, as can be seen 
in Figure 2. The first condition in Eq. Q, $ > 1, is 
equivalent to the requirement that 2?c < given pre- 
viously. From the dispersion relation, Eq. (|3]), we see 
that the resulting destabilization is balanced by the sta- 
bilizing actions of bacterial reproduction and the surface 
tension at large and small wavelength, respectively. The 
unstable modes thus lie within a band qi < q < q2 where 
<?! ^ 9q = v^a/|I?o(po)| and q2 ~ q^ ^ ^|2?e(po)|/K set 
the wavelengths below and above which the stabilizing 
effects of bacteria reproduction and the surface tension 
can compete with the destabilizing effect of the negative 
diffusivity, respectively. For unstable modes to exist, one 
needs qi < q2', restoring prefactors, this yields 2qa < 
which is the second criterion in Eq. Q. This analysis is 
consistent with the view that phase separation is arrested 
by the birth/death dynamics, which stabilizes the long 
wavelength modes (Aq = —a), while the phenomenolog- 
ical tension parameter k primarily fixes the interfacial 
structure of the domains, not their separation. 

We now consider more closely the parameters control- 
ling the transition to pattern formation. For definite- 
ness, we address the specific case used for our simulations, 
-D(p) = £'(0) exp(— Ap). To put Eq. ([T]) in dimensionless 
form, we define rescaled time, space and density as 



t = at: r — 



1/4 p 

r; u = — 
Po 



(5) 



^Pq 



A„ 



q^VM ~ nq" (3) 



The equation of motion now reads 

u = V • [i?e-2*"(l - $u)Vu] + u{l -u)- V^u (6) 

where R = Dt^j^faR and $ = \p{)/2 are the two re- 
maining dimensionless control parameters. Meanwhile 
the conditions Q for pattern formation become 

exp(2$) 



$ > 1; 



R> Rr^2- 



$ - 1 



(7) 
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These relations, combined with the preceding hnear sta- 
bihty analysis, define a phase diagram in the (i?, $) plane 
(Figure [3]) that agrees remarkably well with numerical re- 
sults for systems prepared in a (slightly noisy) uniform 
initial state. 
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Figure 3: Top: Phase diagram in the (R, $) plane. The 
outer region corresponds to stable behavior whereas within 
the curve, patterning occurs. The solid line is the theoretical 
phase boundary - Eq ([7| - which accurately fits the numerics 
(black squares). The blue and red sections correspond to 
continuous and discontinuous transitions respectively. The 
two magenta dots correspond to two 2D simulations which 
show ordered harmonic patterns close to supercriticality and 
amorphous patterns otherwise. Bottom-left: Transition in 
the supercritical regime. The blues lines correspond to the 
theory - Eq ([sjl - whereas the squares comes from simulations 
($ = 1.5; 1.35; 1.2 from top to bottom). Bottom-right: 
Transition in the subcritical regime for $ = 1.06 and $ = 1.7 
(bottom to top). 

Close to the transition, the emergent steady-state pat- 
tern can be studied using an amplitude equation (see 
Supporting Information). Introducing e = {R — Rc)/Rc, 
one gets in ID that for 1.08 < $ < 1.58, the transition is 
supercritical (continuous) and the steady state is given 

by 

M ~ 1 -f A{e) cos{x); 

= , 18(1-'^)^ 

^ ' 34$4 _ 56$3 _ 24$2 _^ 31$ + 19 

which agrees with simulations (Fig 3, bottom left). Out- 
side this range, the transition becomes subcritical (dis- 
continuous. Fig 3, bottom right) and the analytical tools 



available become less reliable [55]. We emphasize again 
that the basic mechanism for patterning presented above 
does not depend on the precise form chosen for v{p). 
Quantitatively however, Eq. (|8]), and the frontier between 
subcriticality and supercriticality, do depend on the de- 
tails of the interplay between the nonlinearity in v{p) 
and the logistic growth term. We leave further analysis 
of such model-specific features to future work. 

While the amplitude equation is more easily devel- 
oped in ID, the stability analysis offered above is valid 
in higher dimensions and it is natural to ask what hap- 
pens in 2D, which is the relevant geometry for Petri dish 
studies with growing bacterial colonies. Fig. 4 shows the 
simulated time evolution of p(r, t) for a system started 
with small random fluctuations around the equilibrium 
density po, with other parameters as in Fig. 1. Per- 
haps not surprisingly, bands are replaced by droplets of 
the high density phase dispersed in a low density back- 
ground at large times. This is the typical steady state 
obtained with a near-uniform starting condition. How- 
ever, the structure and organization of the bacterial drops 
in the steady state depends on the point (i?, <&) chosen in 
the phase diagram. Generally, the closer the system is to 
the supercritical instability curve, the more ordered the 
patterns. For instance we have observed an essentially 
crystalline distribution of bacterial drops, which devel- 
ops defects and eventually becomes amorphous on mov- 
ing further away from the phase boundary (Fig. 3, insets 
to main panel). For particular choices of parameters, our 
model can also admit other steady state patterns. Close 
to the supercritical line, where the phase transition is 
continuous, we can obtain long-lived stripes, whereas for 
fixed large values of R and $ close to the (right) sub- 
critical phase boundary, we have also observed 'inverted 
droplets' with a high density lawn punctuated by low 
density 'holes'. 

In these 2D geometries initialized from a near-uniform 
state, droplets can coalesce in the early stages, while 
at late times the dynamics is governed by evaporation- 
condensation events (see Figs. 4 and 5). However it is 
already apparent from Fig. 4 that coarsening eventually 
stops and the droplets reach rather well-defined steady 
state sizes and centre-to-centre distances. This can be 
quantified by looking at the time evolution of the char- 
acteristic domain size, L(t), which we have computed as 
the inverse of first moment (times 27r) of the structure 
factor 20J. Fig. 5 suggests that L(t) at late times even- 
tually stops increasing and reaches a steady state value. 
(The visible steps in domain size mark discrete evapo- 
ration events involving smaller bacterial droplets; pre- 
sumably L[t) would become smooth for a large enough 
system.) 

These droplet patterns in steady state are very simi- 
lar to those observed for E. coli in a liquid medium or 
S. typhimurium in semi-solid agar (0.24% water-agar in 
Ref. [5]) when starting from a uniform distribution [6]. 
For the E. coli case, interactions are believed to come 
from chemoattractant, emitted by the bacteria them- 
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Figure 4: Numerical results for a 2D simulation with size equal to 16x16, A = 0.3, Do = 1, a = 0.01, n 
Times corresponding to the snapshots are (in simulation units, from left to right): 5, 12, 19.6 and 50.6. 



0.001 and po = 10. 
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Figure 5: Plot of the characteristic domain size, L{t) as a 
function of time for a system in the inhomogeneous phase, 
with initially random density fluctuations around po. Param- 
eters were: a = 0.01, A = 0.27, k = 0.001, while the system 
size was 20 x 20. The solid line corresponds to a single run, 
while the dashed line is an average over 6 runs. The steps in 
the single run curve correspond to evaporation-condensation 
events, highlighted by black squares in the snapshots shown 
in the figure (before and after one of the steps respectively, 
arrows indicate positions on the plot corresponding to the two 
snapshots). 



selves, that is not degraded over time [H H]. The 
chemoattractant distribution should become more and 
more uniform so that these interactions decay to zero 
as time proceeds. In our framework this is analagous 
to decreasing $, which will turn any initially unsta- 
ble state into a homogeneous one, and can thus explain 
that the patterns observed experimentally fade with time 
(whereas in our simulations $ remains constant and the 
pattern are stable indefinitely). E. coli in a semi-solid 
medium also exhibits droplet patterns of high symmetry. 
In our framework, such patterns result from a continuous 
transition, close to the supercritical line. 

The growth of bacterial colonies of S. typhimurium 
starting from a small inoculum of bacterial cells in semi- 
solid agarleads to quite specific (transient but long-lived) 
patterns, with the bacteria accumulating in concentric 
rings that can subsequently fragment into a pattern of 
dots [SlISI- Once again, although these patterns are be- 
lieved to stem from a chemotactic mechanism j6], we 
find they can arise in principle without one, so long as 
our two basic ingredients of density-suppressed motility 
and logistic growth are both present. Indeed, initializing 



our simulations with a single small droplet of high den- 
sity yO, we find that a similarly patterned bacterial colony 
structure develops. First, the bacteria spread radially 
(through a Fisher-like wave), forming an unstructured 
lawn with the highest density at the center. This back- 
ground density increases logistically until the onset of 
instability via our generic phase-separation mechanism; 
with circular symmetry, the instability causes concentric 
rings of high bacterial density to successively develop that 
are very stable in time (Figs. 6c and 6d). The patterns 
observed at later times again depend on the position of 
the parameters in the (i?, $) plane. If we fix a value of 
i?, e.g. 100, larger values of <i> in the unstable region 
lead to rings being very stable. For smaller values of 
on the other hand, effectively corresponding to weaker 
interactions between the bacteria, we observe that rings 
initially form but rapidly destabilize through a secondary 
modulation of the bacterial density along them. This 
eventually breaks the rings into a series of drops. The 
inner rings destabilize first, and the system evolves even- 
tually to the same steady state as found starting from 
a uniform density, composed of drops with well defined 
characteristic size and separation. All this phenomenol- 
ogy is strikingly reminiscent of the dynamics observed 
by Woodward et al [S] for S. typhimirium, where rings 
are stable at large concentrations of potassium succinate 
(a 'stimulant' which promotes pattern formation), but 
break up into drops at smaller ones. Our model shows a 
similar morphological change when decreasing <I>, i.e. the 
strength of the interactions. 

Different views are possible concerning the ability of 
our generic model to reproduce the observed chemotac- 
tic patterns of E. coli and S. typhimurium [6]. One pos- 
sibility is that Eq. ([T|, with the interpretation we have 
given for it, actually does embody the important physics 
of pattern formation in these organisms. Indeed it is 
well accepted that bacteria in the high density concen- 
tric rings are essentially non- motile [27] . The precise 
mechanism leading to this observation is not well under- 
stood [6 , but it is possible that the chemotactic mecha- 
nism mainly acts to switch off motility at high density. If 
so, by focussing solely on this aspect (with a correspond- 
ingly vast reduction in the parameter space from that 
of explicit chemotactic models [3 [6]) our model might 
capture the physics of these chemotactic patterns in a 
highly economical way. Interestingly, our model is es- 
sentially local, whereas chemotaxis in principle mediates 



Figure 6: Dynamics of formation of "chemotactic patterns" in 2D, starting from a single small bacterial droplet in the middle 
of the simulation sample. Top row: Formation of "chemotactic rings" in a system with a = 0.1, A = 0.33, and n = 0.001. The 
simulation box has size 40 x 40. The snapshots correspond to times equal to (in simulation units, from left to right) 10, 50, 
100 and 270. Bottom row: Breakage of "rings" into "dots". The four snapshots correspond to the time evolution of a system 
with a = 0.1, A = 0.26, and k = 0.001. We show a 40 x 40 fraction of the simulation box, with the boundaries far away and 
not affecting the pattern. The snapshots correspond to times equal to (in simulation units, from left to right) 10, 70, 290 and 
1220. 



interactions between bacteria that are nonlocal in both 
space and time. It is not clear whether such nonlocal- 
ity is essential for the chemotactic models in [3 [6] or 
if fast-variables approximations and gradient expansions 
would reduce these models (which invole between 9 and 
12 parameters) into Eq. [1]. In this case, we would still 
have in Eq.[l] a highly economical model for chemotac- 
tic pattern- formation organisms, possibly with a different 
interpretation of I?c and k. 

Alternatively, the success of our local model for these 
chemotactic organisms might be largely coincidental. 
But in that case, such a sparse model should be easily 
falsified, for instance by using the linear stability anal- 
ysis to relate the typical length scale of the patterns 
to microbial parame ters. This length scale is of order 
Stt/qc = 27r-\/|I?o|/a, with |I?e| — D, a typical bacterial 
diffusion coefficient [D - 0(100 ^m^s'^) for E. coli P]). 
Using the previously quoted growth rate a ~ 1 hr"-'^, we 
get a ring separation of '-^ 1 mm, in order-of-magnitude 
agreement with the experimental value [5]. This test 
could perhaps be sharpened usefully by altering the grow- 
ing medium so as to change a. 

More generally, our analysis of Eq. ([T]) shows that the 
main prerequisite for pattern formation, assuming the 
presence the logistic growth term, is negativity of the 
effective diffusion constant P^. For run-and-tumble dy- 
namics, 2?o < was shown to arise for a sufficiently 
strong decay of swim-speed with density; it does so be- 
cause spatial variations in the true diffusivity D{p) cre- 
ate a drift flux pV — —pD'{p)V p/2 which can overcom- 
pensate the true diffusive flux —DWp [T3]. Negative Vc 
could, however, equally arise for any density dependent 
nonequilibrium diffusion process. Indeed, the principle 
of detailed balance, which holds only for equilibrium sys- 
tems, leads to the Einstein relation, that D = kgTfi 



with D a many-body diffusivity and fi the correspond- 
ing mobility. This alone ensures that no drift velocity 
can arise purely from gradients in D. In nonequilbrium 
systems, one should expect generically to find such drift 
velocities, and the run-and-tumble model is merely one 
instance of this. Accordingly one can expect in princi- 
ple to find cases of negative I?e in other microorganisms 
showing distinctly different forms of density-dependent 
self-propulsion. 

To summarize, we have studied the dynamics of a sys- 
tem of reproducing and interacting run-and-tumble bac- 
teria, in the case where interactions lead to a decreasing 
local swim speed with increasing local density. We have 
thereby identified a potentially generic mechanism for 
pattern formation in which an instability towards phase 
separation, caused by the tendency for bacteria to move 
slowly where they are numerous, is arrested by the birth 
and death dynamics of bacterial populations. We have 
shown that these two ingredients alone are enough to 
capture many of the patterns observed experimentally in 
bacterial colonies - including some that have only pre- 
viously been explained using far more complex models 
involving specific chemotactic mechanisms. Indeed, if 
motility decreases steeply enough with density, then a 
spatially homogeneous bacterial population becomes un- 
stable to density fiuctuations leading to the formation of 
bands (ID) or droplets (2D). The length scale of the re- 
sulting pattern is set by a balance between diffusion-drift 
fluxes and the logistic relaxation of the population den- 
sity towards its flxed-point value. Starting instead from a 
small initial droplet of bacteria, we predict the formation 
of concentric rings, each of which may eventually further 
separate into droplets. 

In several well studied systems, such characteristic pat- 
terns are (with good reason) believed to be the direct 
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result of chemotactic behavior [51 [B] . It is therefore re- 
markable that they can also arise purely from the inter- 
play of density-dependent diffusivity and logistic growth, 
without explicit reference to the dynamics (or even the 
presence) of a chcmoattractant. This suggests that sim- 
ilar patterns might arise in organisms having no true 
chemotactic behavior at all. Such patterns could then 
be the result of local chemical signalling without gradi- 
ent detection (quorum-sensing, not chemotaxis) or even 
purely physical interactions (steric hindrance), either of 
which could in principle produce the required dependence 
of motility on density. Last, a motility decreasing with 
density is only one of the many mechanism that could 
lead to 2?o(p) < ^-^^ our analysis would apply equally 
to all such cases. 

The simplest version of our model allows identifica- 
tion of just two dimensionless parameters that control 
the entire pattern-forming process. In both homogeneous 
and centrosymmetric geometries, this gives predictions 
for how the pattern type depends on interaction strength 
which are broadly confirmed by experimental data. This 
suggests that some of the diverse patterns formed by 
colonies of motile bacteria could have a relatively uni- 
versal origin. 
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Appendix A 

We show here how to derive the amplitude equa- 
tion (8). Let us start from the dimensionless equation 
of motion (6) 

ii = V[i?e"2*«(^ _ $u)Vu] + u{l - u) - V^u (Al) 

and recall the two conditions for patterning Eq. (7): 

$ > 1; i?exp(-2$)($ - 1) > 2 (A2) 

To analyze precisely the transition, we derive below the 
steady-state limit of the amplitude equation in ID. By 
inspection one sees that the unperturbed steady-state 
of (All is given by u = 1. To characterize the ampli- 



tude of the perturbation around u = 1, we introduce 
u= 1 + w/^ so that w evolves with 



for e > and small. The dynamics now reads 

Cw~ 2edlw + g{w) (A5) 

where £ = —(1 + 9^)^ is the linear part of the evolution 
operator at the transition, 2ed^w gives an extra linear 
part due to the perturbation (e > 0) and g{w) is the 
non- linear part: 



$ — 1 (j) -^ 

(A3) 

We are interested by the vicinity of the transition where 



w , 



i?e""*($~l) = 2(1-1- e) 



(A4) 



11) r HI 

9H = -^-d. [2(1 + s){{l + — )e- 



1. Amplitude equation 



(A6) 



As usual with the amplitude equation approach, we 
expand w in power series of the perturbation e and study 



Eq. (A5) order by order. As shown below (Eqs (21-23)), 



the correct expansion is 



= + Uie + U2e^^^ + 



(A7) 



Expanding (A6) to the order e'^/^ and substituting 



in (A5| yields order by order: 



CUo = 



(A8) 
(A9) 



-CU, = -2dlUo - ^ - ^29^C/o^AlO) 



4$ 



3$- 1 



Equation ( A8 ) can be easily solved and yields 



(All) 



The amplitude of the perturbation we are trying to derive 
is thus 2\A\. Equation (A9) can also be solved directly: 



Ui = Be" + B*e 



C + De^ 



D*e- 



(A12) 



where B can de determined from higher order equations 
(but does not interest us here), and C and D are given 
by 



C = 



2\Al 



^--(4^--) (A13) 



as can be checked by direct substitution in Eq. ( A9 1 



Equation (AlOl does not always have a solution. In- 



deed, the application of C to any function U2 cannot 
yield a multiple of e", (since £e'^ = and C is linear). 
The r.h.s. however does contain a multiple of e*^ whose 
prefactor must thus vanish. This gives a condition for 
the expansion to provide a proper steady-state solution 
of the problem. Let us summarize the contributions of 
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the different terms to the prefactor of e*^ in the r.h.s of ^fe. One could look for a more general expansion: 



Eq. (AlOl 



3-2$ 
$- 1 



29, 



-2dlUo yields 2A (A14) 
^)UoU, yields 2A|Ap(^l^ - 1()A15) 

13-2$ 
9^ $-1 



w = C/oe" + Uie^ 



(A20) 



4$ 



3$-l 



^a2[/3 yields 4* 



$- 1 



(A16) 



In this case, the expansion of equation ( A5 1 yields two 
power series: ^ RkS"'^ and J2 WkS"'''^^. For the two se- 
ries to give terms that can balance each-other, one needs 
a + I = ka for k > 2 and thus 



The sum of these terms vanishes only if 

y4(^9$2($-l)2-|-2|Ap(34$''-56$^-24$2-F31$-hl9)) = 

(A17) 

and thus either 



1 

k - 1 



(A21) 



0; 



or 



9$2(i - $) 



The candidates for a are thus 1; 1/2; 1/3; .... Note that 
a < 1 implies 2a -I- 1 > 3a. We can therefore stop the 
expansion at 3a and 2a -I- 1 to get the first three terms 



2(34$4 - 56(()3 _ 24$2 + 31$ + igyn the expansion of equation ( |A5 | 



(A18) 

Finally, the first order in the amplitude equation yields 



Let us first try a — 1. The order by order the expan- 
sion yields 



j{x) — 2\A\y/e cos{x — xq) 



(A19) 



where xq is a constant. Note that by construction 
\A\'^ > and a non-zero solution only exists for $ e 
[1.08439,1.59237]. For these values of $, Eq. (tATsl) 



0, 0{e) 



2$, 



$- 1 



(A22) 
0(24)V23) 



and (A19) work very well, as can be checked in figure 



[7j Outside this range the transition becomes subcritical 
and the standard approach does not work anymore. Al- 
ternative treatments have been proposed but are not as 
reliable (see ref [20] for more details). Interestingly, we 
see that the order of the transition and the amplitude of 
the perturbation depend on how non-linear terms in g{w) 
balance the linear growth term —2£d^w in (A5). Since 



Equation ( A22 ) yields C/q = Ae^ +A*e '^"^ but equation 



the former depends on the non-linear relation we 
do not expect equation (AI8) to be generic, as opposed 



to the stability analysis which can be expressed solely in 
terms of I?e(/3o) and its derivative. 



2. What is the correct expansion? 



In ( A7 1 , we expanded w in power series of ^/e, thus 



assuming that the amplitude is an analytic function of 



( A23 1 cannot be solved since there is a non-zero multiple 
oFe*^^ on the r.h.s. {—2d^Uo) which cannot result from 
the application of to any function. Thus a — 1 is not 
an option. 

For a < 1/3, then a -I- 1 > 1 > 3a. There is thus 
no contribution of —2£d^w to the first three orders in 
the expansion of (|A5[ ). In particular, the two first order 
are still given by ( |A8[ ) and (A9), whereas the third order 
is given by (AlO) without the term linear in Uq. This 
means that the contribution ( A14[ ) is not present and 
the prefactor of e'*^^ in the r.h.s. of (AlO I only contains 
multiples of The resolvability condition (A17I is 

thus of the form A\A\'^ f{^) = which implies 0. 
The only expansion which yields a result is thus for a = 
1/2. 
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